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We present a theoretical study on the electron transmission through the AB-BA stacking boundary 
in multilayer graphenes. Using the tight-binding model and the transfer matrix method, we calculate 
the electron transmission probability through the boundary as a function of electron Fermi energy 
in multilayers from bilayer to five-layer. We find that the transmission is strongly suppressed 
particularly near the band touching point, suggesting that the electronic conductivity in general 
multilayer graphenes is significantly interfered by stacking fault. The conductivity suppression by 
stacking fault is the strongest in the bilayer graphene, while it is gradually relaxed as increasing the 
number of layers. At a large carrier density, we observe an even-odd effect where the transmission 
is relatively lower in trilayer and five-layer than in bilayer and four-layer, and this is related to the 
existence of a monolayer-like linear band in odd layers. For bilayer graphene, we also study the 
effect of the perpendicular electric field opening an energy gap, and show that the band deformation 
enhances the electron transmission at a fixed carrier density. 


I. INTRODUCTION 

After the successful isolation of individual graphene 
flake, graphene and its multilayers became one of the 
most intensively studied topics. In multilayer graphene, 
the weak van der Waals interlayer interaction allows do¬ 
main structure consisting of different stacking regions. 
Recent experiments found that bilaye r gr aphene samples 
commonly contain domain boundariesU21 which connect 
two distinct Bernal-stacking structures referred as AB 
and BA (Fig. [I]). In trilayer, similar stacking fault forms 
ABA-ABC stacking boundary which connects the Bernal 
stacked and rhombohedral stacked trilayer graphene®. 
This kind of domain structure is expected to rather com¬ 
mon in graphene multilayers, and it should influence the 
electric conductivity through the scattering of charge car¬ 
riers. The effect of the domain boundary on the elec¬ 
tron transport was studied for graphene bilayei^l, and 
it was found that the transmission through the boundary 
is highly suppressed at the low carrier concentration. 

In this paper, we extend the previous analyses to gen¬ 
eral graphene multilayers from bilayer to five layer and 
investigate the characteristic boundary effect on the elec¬ 
tron transport in increasing the thickness. Specifically 
we consider a multilayer system in Fig. [2] having AB-BA 
boundary on the bottom layer while otherwise stacked in 
Bernal stacking order. We calculate the electron trans¬ 
mission probability through the stacking fault as a func¬ 
tion of the Fermi energy using the tight-binding model 
and the transfer matrix method. The result shows that 
the boundaries generally suppress the electron transmis¬ 
sion in comparision with the perfect Bernal stacked sys¬ 
tem, while the reduction becomes gradually weaker as 
increasing the number of layers. We also notice that the 
transmission shows some odd-even effect at high elec¬ 
tron density larger than 10 13 cm -2 , where the even-layer 
boundaries have better electron transmission than odd- 
layer boundaries. We present the qualitative explanation 
of those characteristic properties of electron transmission 
through domain boundary by considering the local band 


structure inside the boundary region. 

For Bernal stacked bilayer graphene, we can open an 
band gap between conduction band and valence bands 
by applying an external electric field perpendicular to 
the graphene sheet PEU Here we also study the electron 
transmission in bilayer AB-BA boundary under the per¬ 
pendicular electric field, and find that the band deforma¬ 
tion induced by the electric field significantly enhances 
the electron transmission at a fixed carrier density. 


II. FORMULATION 
A. Atomic structure 

Let us consider the AB-BA domain boundaries in bi¬ 
layer graphene as illustrated in Fig. [l] where the uni¬ 
form BA stacking in the left region (x < 0) is continu¬ 
ously deformed in the boundary region (0 < x < W) to 
connect to the uniform AB stacking in the right region 
(x > W). We assume that the system is periodic in y- 
direction. The side view shows the atomic positions on 
a ccz-plane, where filled circles and open circles indicate 
A and B atoms, respectively, while dots are the centers 
of hexagons denoted as C sites (no atoms). The domain 
boundary is created by shifting two graphene layers in 
the right region with respect to each other, while fixing 
the left region. Specifically, we start with the regular 
BA bilayer graphene [Fig[I] (a)], and shift the top layer 
and the bottom layer of the right region by ±a/(2v / 3) 
(where a = 0.246 nm is the graphene’s lattice constant), 
respectively, with the structure in the left region fixed. 
We assume that the atomic positions in the intermediate 
region linearly interpolates between left and right. The 
width of the intermediate region, W is written in terms 
of an the number of unit cells N as W = (N + I) y/3a. 
At N = 20, for instance, we have W ~ 8.7 nm, while the 
realistic boundary width is about lOnnP. Here we focus 
on the boundary parallel to the zigzag direction which 
has the biggest effect on the transmission.P 
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FIG. 1: Atomic structure of (a) regular BA bilayer graphene and (b) BA-AB domain boundary on bilayer graphene. 
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FIG. 2: Side view of N-layer domain boundary, which consists of A — 1 Bernal stacking layers (red layers) and one distorted 
layer at the bottom (blue layer). 


The construction of the domain boundary in general 
iV-layer graphene is similar to the bilayer case, where we 
assume only the bottom layer is expanded by a/ a/ 3 in 
£-axis while other N — 1 layers are fixed to the Bernal 
stacking without distorion as illustrated in Fig. [5] The 
stacking structure can be specified by the sequence of 
sites along a vertical line; here we see that • • • BABA 
stack at x = 0 is continuously deformed to • • • BABC 
stack at x = W. 


B. Tight-binding model 

To model the motion of an electron in the system, we 
use the tight-binding method. The Hamiltonian for tight- 
binding model is written as 


H = -J2 - R j)\ R i) ( R j I + h.c. (1) 


where Ri is the atomic coordinate, |J?i) is the wavefunc- 
tion at site *, and f(-R; — Rj) is the transfer integral 
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between atom i and j. We adopt an approximatiorP^211 


- t(d) = V pp7r (d) 


1 - 


d • e. 


T Vppc t (d) 


d ■ e. 


( 2 ) 

v PP A d ) = V p P K ex P (^^7^) > ( 3 ) 

v ppa (d) = V° pa exp ( 4 ) 


where d = Ri — Rj is the distance between two atoms. e z 
is the unit vector on z axis. V pp7T « —2.7eV is the transfer 
integrals between nearest-neighbor atoms of monolayer 
graphene which are located at distance ao = a/\/3 ss 
0.142nm. V pp<7 ~ 0.48eV is the transfer integral between 
two nearest-vertically aligned atoms, do « 0.335nm is 
the interlayer spacing. The decay length ro of transfer 
integral is chosen at 0.184a. At d > v/3a, the transfer 
integral is very small and can be neglected. 


C. Transmission probability 


Fig.[3]shows the schematic view of Fermi surface in the 
left region and the right region at a fixed energy e. For k y 
falling inside the Fermi surface, we have traveling modes 
with real k x , which can be classified as left-going and 
right-going modes depending on the expectation value of 
velocity in x. Since the system is translationally symmet¬ 
ric along y-direction, the wavenumber k y is conserved, so 
that an electronic state in the left region is only con¬ 
nected to one in the right region on the same horizontal 
line. Using the formulation in the Appendix, we calculate 
the transmission coefficient t pl ,[ky', e) from incident chan¬ 
nel v in the left region to out-going channel fi in the right 
region. The case (a) in Fig. |3]shows a situation where we 
have a single channel per direction i.e, we have the right 
going mode 1 (1') and the left going mode 2 (2') in the 
left (right) region. The electron transmission probability 
is then given by |ti'i| 2 - In (b), we have multiple channels 
per direction; two right going modes 1, 2 (1', 2') and two 
left going modes 3, 4 (3 7 , 4') in the left (right) region. 
We then have multiple transmission probabilities |ti'i| 2 , 
|^ 2 'i| 2 , |U' 2 | 2 , and |f 2 ' 2 | 2 - 

To characterize the total transmission for given bound¬ 
ary structure, we calculate the Landauer’s conductance 


G{e) 


ky PV 



( 5 ) 


where the indexes /r and v run over all the traveling 
modes, and we assumed the system is periodic in y- 
direction with a sufficiently long period L y . The con¬ 
ductance G is naturally proportional to L y , and the con¬ 
ductance per unit width, G/L yi can be used as a quantity 



FIG. 3: Schematic view of Fermi surfaces of the left and right 
regions and traveling modes at a fixed energy and k y . 
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FIG. 4: Energy band structures of bilayer graphene in the 
vicinity of K-point with the interlayer asymmetric potential 
(a) A = 0 and (b) A = 0.1 eV. The dashed lines indicate the 
Fermi levels corresponding to several carrier densities n s . 


which measures the transparency of the domain bound¬ 
ary for the electronic transport. 


III. NUMERICAL RESULT AND DISCUSSION 

A. Bilayer graphene and the effect of asymmetric 
electric potential 

We first consider the AB-BA domain boundary in bi¬ 
layer graphene [Fig. [l]. The width of the boundary region 
is chosen as N = 20, or W = (N + l/2)ai/3 ~ 8.7nm. 
Fig. [I] (a) shows the band structure of uniform AB-bilayer 
graphene, where the energy origin e = 0 is set to the band 
touching point. In the top panels of Fig. [5] (a), we plot 
the transmission probability as a function of wave vector 
k y around the Tv-point at several Fermi energies. Here 
we choose the energy range 0 < Ef < 0.15eV, where we 
have at most a single traveling channel in one-direction at 
each k y , where we have the only transmission probability 
Ill'll 2 . 

We observe that the electron transmission is very small 
near the bottom of conduction band {ef = 0.03eV) where 
the electron waves mostly get reflected at the boundary, 
while the transmission probability sharply rises up in in- 
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(a) A = 0 (b) A = O.leV 



FIG. 5: Electron transmission probability k y )| 2 through AB-BA domain boundary (W = 8.7 nm) on the bilayer graphene 
with the interlayer asymmetric potential (a) A = 0 and (b) A = 0.1 eV. The lower panels show the local Fermi surfaces and 
their corresponding local atomic structures. 


creasing the Fermi levelP^This behavior can also be seen 
in Fig. [6] which plots the total transmission G/L y as a 
function of electron density. In Fig. [bj we also compare 
the total transmission for the AB-BA domain boundary 
(blue) with that for the uniform AB bilayer graphene 
without domain boundary (black). We see that, that ra¬ 
tio of the former to the latter is almost zero near the 
zero energy, while rapidly increases at higher energy lev¬ 
els. This result suggests that the existence of the domain 
boundary dramatically reduces the conductivity of the 
bilayer graphene near the charge neutral point. 

These characteristic features can be roughly explained 
from two different aspects. First, we generally expect 
a bigger transmission in a smoothly changing boundary 
than in the sharply changing boundary. The smooth¬ 
ness / sharpness is characterized by the ratio of the wave 
length of the incident wave A to the width of the bound¬ 
ary region W. In the low Fermi energy (£p = 0.03eV), 
the effective wavenumber measured from K point is 
about k ss 0.04a _1 , and then the wavelength A = 2ir/k ~ 
40 nm is actually larger than the width of boundary re¬ 
gion at W ~ 8.7 nm (N = 20). So this case is regarded as 


a sharp boundary, suggesting that we have a considerable 
reflection. In increasing the Fermi level, the ratio A/IT 
simply decreases, i.e., the boundary becomes relatively 
smoother and gives a larger transmission. 

Another important aspect is the wavenumber matching 
in local Fermi surfaced In Fig. [2j we see that the stacking 
structure gradually changes from AB to BA in increasing 
x, and the local structure at any particular x approxi¬ 
mates the uniform bilayer graphene with a certain inter¬ 
layer translation. Then we can consider the local Fermi 
surface corresponding to the local atomic structure as a 
function of x, which is shown in the lower panels of Fig. 
[5] At £p = 0.03 eV, for instance, the Fermi surface in the 
intermediate region shifts along the k y direction, leaving 
only a small overlap with initial Fermi surface when pro¬ 
jected to k y axis. As a result, for a traveling wave of k y off 
the overlapping region does not have corresponding trav¬ 
eling waves in the intermediate region, then it does not 
reach the other side. Meanwhile, at higher Fermi levels, 
the intermediate Fermi surface better overlaps with the 
initial one on the k y axis, so that the boundary becomes 
more transparent. 
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FIG. 6: Total transmission G/L y as a function of the electron 
density calculated for bilayer graphene of A = 0 with (blue, 
solid) and without (black dotted) AB-BA domain boundary 
(W = 8.7 nm). The red dashed curve is for the domain bound¬ 
ary in presence of the interlayer asymmetric potential A = 0.1 
eV. 


that the transmission probability is much greater in the 
gapped case than in the gapless case. The total transmis¬ 
sion calculation in Fig. [6] also shows that the suppression 
at low electron density observed in the gapless case is 
considerably relaxed in the gapped case, suggesting that 
applying the external electric field definitely increases the 
electron transmission near the charge neutral point. Here 
it should be noted that the total electron transmission 
can be changed purely by deforming the band structure, 
with the electron density fixed. This feature can also be 
explained in terms of the local Fermi surface structure 
shown in the lower panels Fig. [ 5 ] At n s = 1.0 x 10 12 
cm -2 , the Fermi surface spreads in a wider fc-space re¬ 
gion than in gapless case due to the complex structure 
of band bottom. This results in a shorter Fermi wave¬ 
length and also a better overlap with the intermediate 
Fermi surface, and both of them contribute to a larger 
transmission. 


B. Trilayer, four-layer and five-layer graphene 


The correspondence between the transmission and the 
band structure is only approximate, but the agreement 
should become almost perfect when W becomes much 
larger than the typical wavelength. To see this, we cal¬ 
culate the transmission probability at various widths of 
boundary region, W. In Fig. [TJ we show the transmission 
probability in the top panel, and the local band structure 
in the bottom. The blue shading indicates the region of 
k v which has traveling modes throughout the interme¬ 
diate region. For small boundary region W = 8.7 nm 
(IV = 20), we still have a finite transmission probabil¬ 
ity in the region outside of the blue shading region, be¬ 
cause there the wavelength A is actually larger than W 
as already argued, and thus the local Fermi surface is 
not well defined. Increasing W, the transmission out¬ 
side the shaded region vanishes, and at the same time 
the transmission probability in the shaded region rapidly 
rises and reaches almost 1 at W = 85.4nm. This is be¬ 
cause the wave is hardly reflected in a smoothly-changing 
boundary which is much longer than the electron wave 
length, as long as a traveling mode exists throughout the 
intermediate region. 

Next, we apply a perpendicular external electric field 
to the system by adding the asymmetric electrostatic po¬ 
tential +A and —A to the top and bottom layers. As 
shown in Fig. [4] (b), A actually opens up a band gap 
between the touching conduction and valance band, and 
near the band edge, the energy dispersion exhibits a com¬ 
plex feature where the Fermi surface split into several 
different parts due to the trigonal warpingP In Fig. [ 5 ] 
(b), we compare the transmission probabilities in gapped 
case (A = 0) and that in gapless case (A = 0.1 eV) at 
the two different electron densities. At the lower density 
n s = 1.0 x 10 12 cnG 2 , we have two right-going channels 
in the gapped case so that we plot multiple transmission 
probablities |H'i| 2 , |f 2 'i| 2 = |fi' 2 | 2 , and |t 2 ' 2 | 2 - We see 


In more than three-layer stack, the left and right re¬ 
gions have generally different energy bands as shown in 
Fig. [I When we consider a realistic experimental situa¬ 
tion with a single gate electrode underneath the sample, 
we hav e a h omogeneous carrier density over the whole 
system p2124l where the left and right regions and share 
the same Fermi energy, while the origins of band en¬ 
ergies are relatively shifted to achieve the same charge 
density in the both regions. The situation is schemati¬ 
cally illustrated in Fig. [9] for the trilayer case. The shift 
of the band energy is given by the different electrostatic 
potentials between the two regions, which automatically 
arises to satisfy the electrostatic equation in the field- 
effect transistor geometry. In the intermediate region, 
we simply assume the electrostatic potential V ( x ) lin¬ 
early changing in the intermediate region to connect left 
and right as shown in the lower panel of Fig. [9j i.e., 

(0 (x < 0) 

V(x) = { ^ (0 <x<W) , (6) 

I (x> W). 


We calculate the transmission probability for the tight- 
binding Hamiltonian with V(x) is included. 

Figures 10 (a) (b) and (c) plot the total transmission 
G/L y as a function of the electron density for trilayer, 
four-layer and five-layer cases, respectively, where the 
black and blue curves are for with and without domain 
boundary, respectively. Figure 10 (d) shows the ratio 
of the transmission with the boundary to without the 
boundary. The overall behavior in the trilayer boundary 
is pretty much similar to the bilayer case, i.e., the trans¬ 
mission is suppressed near the charge neutral point and 
gradually rises at the high Fermi energy. In adding more 
layer, however, we see that the reduction near n s = 0 
becomes gradually weaker. It is a natural consequence 
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FIG. 7: (Top panels) Electron transmission probability \t\'\(k y )\ 2 in the bilayer domain boundary with various widths W. We 
fix A = 0 and n s =4x 10 12 cm -12 . The blue shading shows the overlap region of the local Fermi sufaces (see text). 


because the wave amplitude per single layer becomes rel¬ 
atively smaller in a larger stack so that the transmission 
is less sensitive to the existence of the domain bound¬ 
ary on the surface layer. When it comes to high density 
region, we observe some odd-even effect in Fig. 10 (d), 
where the boundaries with even number of layers show 
relatively better transmission than in the odd layers. It 
is particularly conspicuous when comparing bilayer and 
trilayer, where the transmission at n = 20 x 10 12 cuU 2 is 
almost 90% in bilayer while it is only 60% in trilayer. 

The even-odd characteristics is closely related to the 
left and right band structures presented in Fig. [8] There 
the horizontal dashed lines indicate the Fermi energies 
corresponding to several carrier densities n s . In trilayer 
at n s < 20 x 10 12 cm -2 , for example, we have two right¬ 
going channels v = 1,2 in the left region, and a single 
right-going channel /i = ¥ in the right region. We notice 
that the odd layer always has one more channel in left re¬ 
gion than in right region, and this is due to the existence 
of the monolayer-like band in Bernal stacked odd layer 
grapheneP^FFig. 11 shows the transmission probability 
\t^{k y )\ 2 from the incident channel v (summed over 
out-going channels fi) in trilayer, four-layer and five-layer 


at n s = 20 x 10 12 cm -2 . In the odd-layers, the incident 
carriers from the linear band (green curve) are not very 
well transmitted, and at the same time, the carriers from 
the other non-linear bands are considerably reflected only 
in the region of k y where a linear band state exists. This 
suggests that the reflection matrix elements between the 
linear band and the other bands are particularly strong, 
and it significantly interferes the total transmission. The 
strong reflection in presence of the linear band is presum¬ 
ably related to the fact that a linear band state has a 
relatively large wave amplitude on the surface layer com¬ 
pared to other bandsPThis fact explains why the trans¬ 
mission at the high carrier density is generally higher in 
even layer cases where the linear band state is absent. 


IV. CONCLUSION 

We studied the electron transmission properties in 
stacking domain boundary in multilayer graphenes from 
bilayer to five-layer. We find the boundary significantly 
reduces the electron transmission at low Fermi energies, 
while the reduction becomes gradually weaker as increas- 
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FIG. 8: Energy band structures in the vicinity of K-point 
of the left region (left panel) and right region (right panel) 
in trilayer (top), four-layer(middle) and five layer (bottom) 
domain boundary. The dashed lines indicate the Fermi levels 
corresponding to several specific carrier densities n 3 . 


ing the number of layers. We also find an odd-even effect 
at high electron density, where the even-layer boundaries 
have better electron transmission than odd-layer bound¬ 
aries. For bilayer, we also demonstrated the conductance 
suppression at low energies is significantly relaxed by ap¬ 
plying a perpendicular external electric held. We found 
that the transmission characteristics is closely related 
the local electronic structure inside the boundary region. 
The agreement between the transmisssion and the over¬ 
lap of the local band structure is generally better for a 
wider boundary, and it becomes almost prefect when the 
boundary width is much greater than the typical Fermi 
wavelength. 

The study of the electron transmission through a sin¬ 
gle domain boundary provides a very simple model to 
roughly understand the transmission properties of macro- 


ABA ABC 



FIG. 9: Schematics of the electrostatic potential V(x) to 
achieve the same electron density in ABA stack and ABC 
stack regions in trilayer graphene. 


scopic multilayer graphenes and bulk graphite, which are 
expected to have a number stacking faults and domains. 
Our result indicates the electron transmission in a gen¬ 
eral graphitic system is significantly supressed in presence 
of domain boundaries. The present result also suggests 
that a boundary between two different stacking struc¬ 
tures may provide a mechanism to control the electron 
current on multilayer graphene. In particular, the fact 
that the system being almost insulator near bottom of 
conduction band and quickly becoming transparent when 
increasing Fermi levels gives a possibility to be applied 
to switching devices. 


Appendix: Transfer matrix method 


In this section, we present the formula to calculate the 
transmission probability in the domain boundary.^ We 
first consider the semi-infinite left and right regions sep¬ 
arately, and specify the traveling modes in each region. 
Then we consider the intermediate region and calculate 
the transfer matrix connecting the traveling modes in the 
left and right. 

In uniform AB-stacked bilayer graphene, for example, 
we define a unit cell as in Fig. [l2j which consists of eight 
atoms. The Schrodinger equation reads 

(e - n 0 )Cj = —T^Cj-i - TC j+1 , (A.l) 

where C :) . Cj_ i, C 7+1 are 8-component vectors consist¬ 
ing of 8 wavefunctions of 8 atoms in each unit cell, and 
w 0 , rt, T are 8x8 matrices. Wo consists of transfer in¬ 
tegrals between atoms inside each unit cells, while T', T 
consists of transfer integrals between atoms in j-th cell 
and atoms in j — 1-th, j + 1-th cell respectively. The 
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FIG. 10: Total transmission G/L y as a function of the electron density in (a) trilayer, (b) four-layer and (c) fivelayer, where 
the black dotted and blue sold curves are for with and without domain boundary (W = 8.7 nm), respectively, (d) Ratio of G 
with the boundary to that without the boundary from bilayer to five-layer. 


equation is transformed as 


(A.2) 

which gives 16 eigenvalues of A. These inodes are clas¬ 
sified as traveling inodes if |A| = 1, and as evanescent 
inodes if |A| ^ 1. For traveling modes, we can define the 
real wave number k x by A = exp (ik x \/?>a), and it is even 
classified to left-going or right-going modes depending on 
the expectation value of the velocity along cc-direction. 
For evanescent modes, we categorize |A| > 1 and |A| < 1 
as left-going and right-going modes, respectively, for the 
sake of convenience. The 16 eigenmodes are always com¬ 
posed of 8 right-going and 8 left-going modes. Let A±, 
Xf , ..., Ag be the eigenvalues of the right and left going 
modes, respectively. 

The wavefunction is written as Cj = A J u with 8- 
component eigenvector u. We define uf as the eigen¬ 
vectors corresponding to Xf. We then define the 8x8 
matrices 


A ± = 




(A.4) 


Now any left and right-going waves can be written as 
the superposition of the waves functions uj or u+. For 
example, wave function at site 0 is written as. 

8 

Co = J2 u f a ' = u±a± > ( A - 5 ) 

i 

where 


(A.6) 


Then the wave amplitude at j -th site becomes 



C f = E(A ^nfaf = £/ ± (A ± y« ± , (A.7) 

i 

leading to the relation 


(A.3) Cf + i = [^(A ± ) J '(I7 ± )- 1 ]C7± = F ± Cf, 




(A.8) 
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FIG. 11: Transmission probability \t^v(k y )\ 2 through the 
domain boundary (W = 8.7 nm) from the incident channel 
v in (a) trilayer, (b) four-layer and (c) five-layer at n s = 
20 x 10 12 cm -2 . The channel index v = 1,2,• • ■ corresponds 
to the numbers in Fig. [8] 



FIG. 12: Unit cell of AB-bilayer graphene. 
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FIG. 13: Schematic view of the AB-BA domain boundary. 


which is similar to the Bloch theorem. The wavefunction 
at site j can be generally written as summation of left¬ 
going wave and right-going wave. 

C j = C+ + CJ (A.9) 


By doing the similar process for BA bilayer graphene 
(the right region), we are also able to get the eigenvec¬ 
tors and eigenvalues. We define U^, A*, for the left 
region (AB-stack), and U^, A^, for the right region 
(BA-stack). 

Figure |~Pj| schematically illustrates the Hamiltonian of 
the AB-BA domain boundary. The tight-binding equa¬ 
tions read 


(e - Ho)Co = -ToCinter - T^C-i, 

(e — Winter Winter = — Tg Co — FjV+lCjv+l , 

(e — 'Hn+i)Cn+1 = — Fjv +1 C'i n ter — TrCn+2- 

(A.10) 

Here Co, C_i, Cy+i, Cjv +2 are 8-component vectors, 
Ci n ter is n-component vector, with n being the number 
of atoms in the intermediate region. Ho, Hn+i are 8x8 
matrices, H ln ter is a nxn matrix, and To, 7/v+i are n x 8 
and 8 x n matrices, respectively. 

In the following, we consider a situation where the in¬ 
cident wave comes from the left region, and transmits to 
the right or reflects back to the left. From (A.8), we have 


C-i = (F+)- 1 C 0 + + (FJ ) _1 Cg" 

= (F+)^C+ + (FT)- 1 (Co - C+), (A.11) 

Cjv +2 = FftCx +1 = Fr Cjv+i- 


Note that there are only right-going waves in the right 
region. Then (A.101 becomes 


f (e-Ho + TliF+y^Co 

I = -T 0 C in ter - ^[(itf)"*) ~ (F£)~ 1 }C+ , 

| (e — 'Hi„ter)Ci n ter = — Tq Cq — TjV+1 CjV+1, 
l (£ T.N+1 + TrF^)Cn + i = —Tjj +1 Ci n te r . 

(A.12) 

which leads to 


/Co 

I Cj n ter 
\Cy + i 


1 

E-H 


-Tt[(F+)-i) - ( f -)- 1 ; 


C+ 


0 

0 


(A.13) 
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The Hamiltonian is (n + 16) x (n + 16) matrix 



fn 0 - 

- T 1 2 ( F +)- 1 

-T 0 

0 \ 

a = 


~ T o 

'Winter 

—Tn+i 


{ 

0 

i 

AS 

55 — 

+ 

Rn+i — TrF+ J 





(A.14) 

We define the Green function 




1 

(G 1 ■■■ \ 



G 

1 

Us ••• J 

, (A.15) 

~ E-n 


where G\ and G 2 are 8x8 matrices. 

The wavefunction at site j = 0 can be written as 

Co = C++C 0 - = -C 2 Tt[(F+)- 1 -(F-)- 1 ]C 0 +. (A.16) 


Then we have 

Co = {- G 2 Tl[(F+r 1 - - 1}C+. (A.17) 

Similarly, the right-going waves at site N + 1 is writte 
nas 

C+ +1 = C N+1 = -G 1 Tl[(F+)- 1 - (F i -)- 1 ]C 0 + . 

(A- 18 ) 

To calculate the transmission and reflection probability, 
we first write the left and right-going waves at site 0, N 


+ 1 in terms of left-going waves and right-going waves 

C+ = U+a t , (A.19) 

C~ = U£a r , (A.20) 

Cft+1 = Ur a t- (A.21) 

The transmission and reflection matrices are defined by 
a t = Tcxi and a r = Rati, and they are given by 

T = (U+r 1 (-Gz.TiKF+y 1 - (F£ ) -1 ]) U+ , (A. 22 ) 

R = (UE )- 1 (-Gi.TlKF+r 1 - (F £) _1 ] - l) U+. 

(A.23) 

The transmission coefficient for incident wave v with ve¬ 
locity v v and out-going channel /.t with velocity v lt is writ¬ 
ten as 

/ j \ 1 ^ 2 3 4 5 6 * 8 9 10 11 12 

(A.24) 

and the reflection coefficient for incoming channel v with 
velocity v v , and out-going channel /.i with velocity i> M as 



These formula can also be used to derive the transmission 
probability for trilayer case. 
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